Motion of charged particles on the Reissner-Nordstrom (Anti)-de Sitter black holes 



Marco 01ivare£| 
Departamento de Fisica, Facultad de Ciencia, 
Universidad de Santiago de Chile, 
Casilla 307, Santiago 2, Chile 

Joel Saavedr£0 

Instituto de Fisica, Pontificia Universidad de Catolica de Valparaiso, 
Av. Universidad 330, Curauma, Valparaiso , Chile 

J.R. Villanuev£@ and Carlos LeivaH 
Departamento de Fisica, Facultad de Ciencias, 
Universidad de Tarapacd, 
Av. General Veldsguez 1775, Arica, Chile 
(Dated: February 8, 2011) 
O ■ 

, In this work we address the study of movement of charged particles in the background of charged 

black holes with non-trivial asymptotic behavior. We compute the exact trajectories for massive- 
charged particles in terms of elliptic Jacobi functions. Finally we obtain a detailed description of 
orbits for Reissner-Nordstrom (Anti) de Sitter black holes in terms of charge, mass and energy of 
the particles. 
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I. INTRODUCTION 



(N 

J> ' Motion of particles on black holes is one interesting phenomena of classical gravity. At this respect, there are several 
OO studies about geodesic motion in the vicinity of black holes. For instance, the geodesic structure of the Schwarzschild 
' (S), Reissner-Nordstrom (RN) and Kerr (K) black holes were studied in detail by Chandrasekhar There, the 
author studied motion of the test particles using the Lagrangian and the Newmann-Penrose formalism. If we think 
this problem in a modern point of view, it is of high interest to study the motion of particles in the vicinity of black 
holes with asymptotic behaviors others than the Schwarzschild case (flat case) . Then, the so called S AdS or SdS black 
holes have been objects of highly consideration, due to the AdS-CFT conjecture, where there are powerful tools that 
relate gravitational theories with an asymptotically AdS behavior with a conformal field theories of one less dimension 
0. Besides, this kind of black holes provides a theoretical laboratory for understanding a lot of important points in 
black holes physics and gravity theories and their extensions. It is possible to generalize the asymptotic flat spacetimes 
\ by including a non-zero cosmological constant term. For instance, the Kotler solution [3|is obtained by imposing 
' R^u = A(7^y in the field equations. Some important implications of those can be found in [^. Different focuses in 
■ the research can be developed depending on the value of the cosmological constant: the Schwarzschild de Sitter case 
(SdS) (which is obtained by considering a positive value of the cosmological constant A > 0), or the Schwarzschild 
anti-de Sitter case (SAdS) (which is obtained by considering a negative value of the cosmological constant, consider, 
A < 0). There are some previous works about the geodesic motion in the SdS case [3-^. On the other hand, when 
a negative cosmological constant, A < 0, is taking in to account i.e. Schwarzschild Anti de Sitter case (SAdS), the 
studies are more complicated. However different aspects of the geodesic structure of this spacetime can be found in 
9|-[lll|. The motion of (neutral) particles in RN black hole with non-zero cosmological constant has been studied in 
Furthermore, in Ref. |13l4l5| . authors presented analytical solutions of the geodesic equation of massive test 
particles in higher dimensional Schwarzschild, Schwarzschild (anti) de Sitter ,S(A)dS, Reissner Nordstrom, RN, and 
Reissner Nordstrom (anti) de Sitter, RN(A)dS spacetimes and they obtained complete solutions and a classification 
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of the possible orbits in these geometries in term of Weierstrass functions. Also, the ecuatorial circular motion in 
Kerr- de Sitter spacetime is studied in [l6j . 

In this article we are interested in the study of the motion of charged particles on the Reissner-Nordstrom (Anti) 
de Sitter black hole. In doing so, we start considering the Hamilton- Jacobi formalism in order to write equations of 
motion and then wc reduce our problem to quadratures. Then we solve quadratures equation and we found all the 
possible orbits for the geometries under consideration in terms of the elliptic Jacobi functions. 



II. MOTION OF CHARGED PARTICLES IN THE VICINITY OF BLACK HOLE 



We are interested in the study of the motion of massive-charged particles near to an static, spherically symmetric 
and charged black hole when a non-zero cosmological constant is taken into account. The A > case is called the 
Reissner-Nordstrom de Sitter spacetime (RNdS), while the A < case is call it the Reissner-Nordstrom Anti de Sitter 
spacetime (RNAdS). 

In terms of the usual Schwarzschild coordinates (t, r, 9, (/)), the metric is written as 



ds^ = -fir)dt^ + -— + r'^(de'^ + sin^ed<f), (1) 



where /(r) is the lapsus function given by 



and the coordinates satisfy the relations: — oo < t < oo, r>0, 0<6'<7ryO<(/)< 27r. 

We would like to consider the motion of test particles with mass m and charge q in the framework of general 
relativity. In order to obtain the equation of motion, we apply the Hamilton- Jacobi formalism. In this sense, the 
Hamilton- Jacobi equation for the geometry described by the metric g^i, is 

where represents the vector potential components associated with the charge of the black hole (because, we are 
considering charged static black holes where the only non-vanishing component of the vector potential is the temporal 
A-o — 5* corresponds to the characteristic Hamilton function. Considering our metric this equation can be written 
as 



1 fdS qQV U^SV 1 fdSV 1 fdSV 2 n 

in order to solve this equation we use the following ansatz 

S = -Et + Si{r) + S2 (0) + J<j), (5) 

where E and J are identified as the energy and angular momentum of the particle. Using this ansatz Eq.Q reads as 
follow 

E+^)\r^f{r)(^)\(^)\,Pcsc^e + r^m^^O. (6) 



/(r) \ r J \ dr J \ d9 

Using the standard procedure we recognize the following constant 

i^=(^)%jW^, (7) 

without lack of generality we consider that the motion is developed in the invariant plane = ^ and in this case L is 
equal to angular momentum J. Then, we obtain the equation of motion 

-E+^]\ r^f{r) (^X + L' + r^m' = 0, (8) 



/(r) \ r J \ dr 
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and thus, we find formal solutions for the radial component of the action 



where e = ±1. Now, using the Hamilton- Jacobi method, we simplify our study to the following quadrature problem 

(10) 

\ ' ~ / 

where we can obtain the radial velocity to respect the coordinate time 



1/2 

dr f qQ ^ ' ' 



r 



Now, the condition of turning point (^)r=rt = allows us to define an effective potential. In fact, considering that 

^^'^^ ^0,Vr, (12) 



we can factorize the term under the square root as: 



E- — ]- m (^' + ^]=iE- V-){E - V+), (13) 



where we can recognize the effective potential for the particle with mass m and electric charge q as 

V± {r) = ^± ^)\fn^~^, (14) 

therefore, ea.([TT|) can be written as 

| = ±7^^v4^^Tqo]W^T^. (15) 



Since the charge of the test particle is much smaller than the mass [q ^ M), eventually V- <Q and for this reason we 
choose the positive branch of the effective potential: V^f f ^ V+ = V . The behavior of effective potential is shown in 
FigUfor RNdS and RNAdS cases with different values of the angular momentum, but kepping m, q fixed. It's worth 
noting that the particle still has electric potential energy at the horizon. 

Finally, we can classify different kinds of motion for massive-charged particles through the values of the L in the 
following way 

• Motion of charged particles with angular momentum 

• Motion of charged particles with vanishing angular momentum 

In the following sections we study in detail the orbits before mentioned. 



A. Circular Orbits 



The orbits can be classified by their values of energy and angular momentum. In order to have at least a stationary 
system the effective potential V{r) have to exhibit extremes for fixed values of radial coordinate, r = r^^ 



(16) 




FIG. 1: Effective potential for charged particles on the geometry of RN for different values of angular momentum, L (left 
panel SAdS case and right panel SdS case). Both of them were plotted considering Q — 0.85M, m — 0.2M, q — O.lSAf and 
e = 400M. 



For simplicity, we rewrite the effective potential as 

V{r) = hir) + Vf{r)j{r), (17) 

where /(r) is the lapsus function, and 



Therefore, using eq. (IT7)) into eq. yields 



h'{r)^f{r)j{r) + i(/'(r) j(r) + /(r) /(r)) 



(18) 



Notice that this equation leads to a polynomial of twelfth grade, so, their solution is restricted to the numeric plane. 
However, it is possible to determine the periods of revolution of the circular orbits, both stable and unstable, with 
respect to the proper time, r, and coordinate time, i, in the following way: the condition (|18p allows to obtain the 
angular momentum for the stable and unstable circle orbits, Lx = Lg and Lx ~ Lc (at Tx = and Tx — Tc in FIG. 
[21 respectively; and rx = and Tx ^ ri,Ri in FIG. [71 respectively). In our case, it is given by 



a + bLl^[c + dLl]^, (19) 

where 



a = m firx)hirx) ; ; c= ; ^=^^5 (20) 

Thus, the real solution of the quadratic equation for Lx give us the angular momentum of the circular orbit. 
Explicitly, the constants of motion, Lx and Ex , for the circular orbits are given by 



Lx — 



'b - 2cd - yPT4d(d^~6^ 



-,1/2 



2d2 



(21) 



and 
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m 



b-2cd- + 4d(da - he) 



(22) 



Therefore, the proper period of the circular orbit (Tt — ~^^) is 



b-2cd- ^\P- + 4d(da - he) 
2d2 



-1/2 



(23) 



and the coordinate period (Tf = l'^s'(j ) ~ hir^)) — TrJ jjy^) is 



b-2cd- y/b^ + Ad{da - be) 
2d?rl 



1/2 



(24) 



Notice that, if neutral particles are taken into account (i. e. g = and h{r,j.) = 0), the proper and coordinate 
period are given by 



Tr = 



2^r, / rl ~ 3Mr, + 2Q^ 
~y Mr,- Q^- Art/3' 



(25) 



and 



Tr 



Wl - 2Mr, + Q^- krtm ^'^ 
rl - SMr, + 2Q2 



(26) 



respectively. Recently, the authors in [17] have studied other aspect of the circular motion of the neutral particles 
in the RN spacetime such as the stability of the orbits. 

III. CHARGED PARTICLES ON THE GEOMETRY OF REISSNER-NORDSTROM ANTI-DE SITTER 

BLACK HOLE 

First at all, let's focus in the motion of charged particles on the geometry of Reissner-Nordstrom Anti-de Sitter 
black hole. In this case charged particles are affected by an effective potential described in right panel of Fig. 1. 
Considering motion with L 7^ we obtain the equation 



f^ = [E-V^ir)][E-V{r)]^~{^y^, 



(27) 



where 



Pair) = J2Pi^^' 



(28) 



and the coefficients are 



Po 



5—, Pi 



-,2 



P3 = -2f M 



qQE\ 



, P5=0, P6 = 1- 
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The solutions (different and reals) of this sixth degree polynomial, correspond to the physical distance that charac- 
terize the motion for the so called periastron, apastron and circular orbit radii. Physical reality of orbits depends on 
the values of constants E, M, Q, q and L. We can distinguish in our description two sets of values for the constants in 
order to classify possible motions. One of them allows the existence of planetary orbits, in this case the polynomial 
has six real distances. Besides, this set has a second class trajectory that represents free fall to the event horizon. 
The second set of fixed values for the constants correspond to critical orbits, for example, here we have the unstable 
circle orbit. Finally, we discuss five trajectories that have physical meaning for the geometry under consideration (see 
FIG. [5] for a better visualization of the orbits.) 

• Planetary Orbit: In this case the orbit corresponds to a bounded trajectory that exhibit oscillation between 
two extremal distances: the periastron and the apastron, {rp and in FIG.[2l respectively). For simplicity, we 
shall consider the case of a double-degenerate real negative root, ap, such that the polynomial Pq can be writen 
as Peir) — {r — ap)'^ Piir). This case has also a particular solution called Reissner- Nordstrom limit. 

• Second Kind Trajectory: This trajectory is computed with the same parameter than the planetary orbits. It 
corresponds to a trajectory that starts at rest from a finite distance, rp. This kind of motion represents the fall 
to the event horizon, and it is considered that has a turning point inside the Cauchy horizon, pp. 

• Critical Trajectory: There are trajectories of the first and second kind with angular momentum Lc, given by 
(??), and energy Ec given taken r = Tc in (jl3p and equating it to zero. The first one starts at rest at a finite 
distance outside from the unstable circular orbit, ri, and then it approximates asymptotically to it. The second 
kind approximates to the unstable orbit from the inside of it. 

• Radial Trajectory: These trajectories have null angular momentum and physically they describe radial fall from 
rest to the event horizon. 

• Circle Orbits: For some fixed values of the constants it is possible to find solutions of the equation of motions 
that represent stable and unstable circle orbits. This case was discussed in the previous section and it is possible 
to find the connection with the Lyapunov exponent for the unstable orbit with the quasinormal modes [cardoso, 
work in progress] . 




FIG. 2: Effective potential for charged particles on the geometry of RNAdS for M = 1, g = 0.85M, £ = 400M, L = 0.33M^ 
q = 0.18M and m = 0.2M. 



A. Planetary Orbit 

In order to characterize this motion we rewrite the polynomial P4:{r) in term of his roots 

P4{r) = (r - ps){r - rs){r - rp){rA - r). 



(29) 
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where ap < < ps < rs < rp < ta, and we identify periastron and apastron distances as rp and ta respectively. 
These trajectories are defined inside the following limits rp < r < va, and the equation of motion has the following 
quadrature 



where Lp is the angular momentum of the test particle in the planetary orbit. The corresponding integral is given by 

Jta [r - (Tp)y/P4(r) 

now introducing constants ap = ^^^^_l'^f^^^_^^y Pp = (^^3^ y 7p = " (rjlt-)(pi-a^) ^ possible to get the 
solution given by a Jacobi elliptic integrals of first and third kinds 

^pM-^r+'^r, (32) 
where 4>p = apf3pF{ipp; Kp) and (ftjj = ap7pn('i/;p; Kp, np), with the elliptic parameters given by 



, . (rp- Ps){rA-r) 

V'p = arcsmW- (33) 

V [rA - rp)(r - ps) 



{ta - rp){rs - ps) 

fip = \ 7 77 r, (-34) 

y (rA - rs)(rp - ps) 

{rp - rA){aF - ps) ^o^A 

np = 7 77 r, (35) 

(rp - ps)\<JF - rA) 

similar results can be found in [isl - fisj . where authors presented analytical solutions for the geodesic equation 
of massive test particles in higher dimensional Schwarzschild, Schwarzschild(anti)de Sitter, ReissnerNordstrom and 
ReissnerNordstrom(anti)de Sitter spacetimes. In FIGl3]we show orbits for an specific value of the angular momentum 
and present the case of an elliptic orbit that precesses between the periastron and apastron. 

At this point, we can calculate the precession of the perihelion as 

20(rp) = 27r + Ayj, (36) 
and considering the exact solution (15^ we obtain 

= 2 [^[f ) (rp) + cp'^P (rp)] - 27r, (37) 

One particular solution is obtained, when we perform a fine tuning in the physical distance under consideration. 
This case corresponds to the called Reissner- Nordstrom limit, that we are going to discuss in the next subsection. 

1. Reissner-Nordstrom Limit 

One approximated solution of the first order represents the limit case of the motion on Reissner-Nordstrom black 
hole. This limit corresponds to the solution when £ — >■ cx) and therefore we must consider the potential 



in which case we obtain two kinds of orbits: 



• confined orbits: they are obtained when the relation < m is satisfied. 
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• non-confined orbits: they are obtained when the relation E > m is satisfied, 
i).- Confined orbits : defining = (m^ — the motion equation (|30p can be written as 

fdrV ^P^jr) 

where, 

4 
3=0 

and the coefficients are 

In terms of the roots of the polynomial,we can write 

P4 = i^A -r){r- r'p){r - r'p){r - p'p). 

Then, the solution can be written as 
where 



a 



p 



Wr'A-r'p){r'p-p'py 



i^'p 



ir'p-p'p)ir'A-r) 
(r'A~r'p){r - p'p) 



, ^ i {r'^~r'p){r'p~p'p) 
Ur'^-r'p){r'p-p'p) 

Defining the constants -00 = -^- and 5 = ['^■^ J't \ , we can write the solution in terms of ehiptic Jacobi sine 

In our case, the precession of the perihelion is given by 

= 2[a'pis:(K'p) -tt], (39) 
where K{k'p) is the complete elliptic Jacobi integral of the first kind. 

B. Second Kind Trajectory 

This trajectory corresponds to the one where the test particle starts at rest from a finite distance bigger than the 
event horizon and then falls to it. Its motion is developed in the following region ps < r < rs- In this case the 
corresponding angular quadrature is given by 



FIG. 3: Polar plot for the planetary orbit of the test particle in the Reissner-Nordstrom Limit. Clearly, there is a precession of 
the perihelion, which is given explicitly by ea. (l39p . 



Jrs [r - aF)\/P[r) 

now we have a fourth degree polynomial described by 

P{r) = {r ~ ps){rs - r){rp - r){rA - r). (41) 

Using a similar procedure than the used for the planetary orbit, we first define the following constants as — 
—, — ap, — — 3 — and 7s = -, — _''-^r/^ _ — r. The solution can be written in terms of Jacobi elliptic 

functions as follows 



Ur) ^ + 4>'i\ (42) 



where (jl-p = as^sF^ips; Kg) and = asjsTl{ips; Ks, ?^s), with 



, . {rp ~ Ps){rs ~ r) 

ips = arcsm J -, (43 

V yrs - Ps)\rp ~ r) 



/ {ta - rp){rs - ps) ^ 

Ks = \ 7 77 r = Kp, (44) 

y [rA - rs){rp - ps) 

Us ^ (;^S-ps)i^P-rp) ^^^^ 
(rp - ps)(crF - rs) 

This kind of trajectory is shown in Fig. (4) where we show the analytic continuation to the inner space to the event 
horizon. In sum this trajectories are doomed to cross the event horizon. 



C. critical trajectory 

There are trajectories of the first and second kind. The first one starts at rest at a finite distance outside from the 
unstable circular orbit and then it approximates asymptotically to it. The second kind approximates to the unstable 
orbit from inner distance. 
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1. First class critical trajectory 



Critical trajectory of the first kind corresponds to the motion of particles to asymptotically tend to a circle orbit 
from a great distance compared to the radius of this orbit. The region of this motion is defined by the limits < r < ri 
and the respective polynomial becomes 



Fc{r) = (r - r6)(r - r5)(r - r4)(r - rcf{ri - r). 
then we can obtain the integral for the orbit as follows 

Defining the following constants ai 



(46) 



(47) 



1 {r - rc)y/{ri - r){r - r4)(r - r5)(r - re) 

^^^^ /3i = ^^^i^^) and 71 = (^i-r^Hre-r,) ^^"^ obtain a solution 



for polar angle 



{ri-rs)(r4-re) ' 



4i)(r)=0«(r)+0W(r), 



(48) 



where the angles 0^'' (r) and cfi'^' (r) are given in terms of the elliptic Jacobi integrals of the first and third kind, 

(l)p\r) = ai/3iF(i/'c; Kc), 0n''('") = ai7in(i/'c; ric), 



respectively, where their parameters are 



ipl.^\r) = arcsin t 



(r4 - re){ri - r) 
[ri - ri){r - rg) ' 



(49) 



(50) 



and 





- r4){r5 


-re) 


jin 


- r5){r4 


-re) 


(re- 


' rc){r4 - 


ri) 


in- 


- rc)ir4 - 


re)' 



(51) 



(52) 



In FIG. 13] we plot the polar angle, 4'i^\ as a function of the radial coordinate, r. We can see that 4>\^' dominates the 
behavior of ipi^^ along the trajectory, specially when r — >■ Tc, in which case (pp^^ takes a finite value. This trajectory is 



(1) 





FIG. 4: Plot of the polar angle, (j>i^\ as a function of the radial coordinate, r, with M ^ 1, Q = 0.85M, ^ = lOM, L = 9M^ 
q = 0.18M and m = 0.2M 



characterized by an asymptotic tendency to the unstable circle orbit from distance bigger than the circle orbit radius 
i.e. in a infinite time we can found the particle in such orbit. 
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2. Second class critical trajectory 



In this case particles asymptotically tend to the circle orbit from a distance less than the critical and the motion 
is realized in the following range < < r < r^- Then, we start considering ai (as was defined in previous case), the 



constant B2 — 1 — ^ — r and 79 = 7 — ^^\r'^ — r, therefore our solution is written as 



4)(r, 

0(2) = ai [^2i^(V'c;Kc) +72n(V'c;Kc,nc)] . 



where 



■0(2) = arcsin , 



(ri - r4){r ~ rs)' 



(53) 



(54) 



(ri - 7-4) (rs - rg) 
(ri - rr,){r4 - rg)' 



and 



^(2) ^ {^5 -rc){ri ~r4) 
{ri - rc)iri - r^) ' 



(55) 



(56) 



As in the last case the motion asymptotically tends to the unstable circular orbit. For example in FIG. [5] we show the 
possible motion starting from a distance less than the one corresponding to the unstable circular orbit and tending 
to it. 





FIG. 5: Second class critical trajectory 



IV. RADIAL TRAJECTORIES 



Radial motion corresponds to a trajectory with vanished angular momentum. First at all, we are considering the 
case of background of RNAdS, where Q < M. Let us start with the relation 



1 



/(r) F{r) ' 



(57) 



where the fourth degree polynomial is given by F{r) = + £2(7-2 — 2Mr + Q^) = (r — r+)(r — r^)j{r). Now the 
dynamical behavior is governed by 



dr 



dr 



qQ 



= S-— -mMl + + 



2M Q2 



£2 



(58) 
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then, the trajectory equation can be written in the following form 



(^) =^^> (59) 



where we define the fourth degree polynomial G{r) = ^j^oCtjr^ , with the coefficients ao = —Q £ {1 — q /m ), 
ai = 2£^{M - qQE/vn?), a2 = ^{E"^ jm^ - 1), as = and a4 = -1. The proper time quadrature is 



TO 



r(r) =±- / 1. dr, (60) 



this motion corresponds to the radial fall starting at rest from the distance i?o- Then, considering negative sign, 
the radial function in this case is written as G(r) — {Rq — r)(r — po){r — oi){r — a^), whose four real roots are 
Rq > r > pq > Q > (Ji > a2- Then we can write our solution as follows 

r(r) = A;on(*; k, u) + fc/F(*; k), (61) 

where 

ko = ± ^ , (62) 

kj = ^ , '"^ (63) 
m y/{Ro - fTi)(po - 0-2) 



,T, • {po - (J2){Ro - r) 

* = arcsm . / — , (64) 

V (^0 - Po)[r - (72) 



and 



' (■R0-P0)(gl -(72) 

{Ro-a,){po-a2y ^ ' 



V = , (66) 

PO - cr2 



are constants. On the other hand, the corresponding coordinate time quadrature is given by 

in order to see physical contents, it is necessary to separate in partial fractions the term under the integral. 



Er^ - qQr^ A B C + Dr 

= ^ ^ (6^) 

/< (r) r — r+ r — r- j{r) 

where function j{r) — + ar + b'^, and the constant are given by a = r_ + r+ y 6^ = + — r_r+. One important 
point is the physical behavior of the coordinate time near the event horizon. There, the leading term is of the first 

order, where the constant is A = . Therefore we can write the general solution in this case as 

t{r) = AotB(r) + Aiii(r) + A2t2{r) + A^^t^ir). (69) 
There, we can see that near horizon the coordinate term is given by the leading term 

t{r) « Aotoir), (70) 
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where the constant is 



2A 



m (i?o - r+){r+ - (J2)\/ {Ro - <Ji){po - (^2) ' 
the exact solution for this term at the event horizons is 

toir) = {Ro - fT2)n(^'; K, n) - (i?o - k), 



(71) 



(72) 



where this solution is defined with the same set of parameters than the proper time solution, but just we are added 
the constant 



r+ - Rq' 



(73) 



FIG. [SI shows the behavior for the proper and coordinate time. In simples terms, these solutions are very similar to 
Schwarzschild case where in the proper time framework particles can cross the event horizon and for the coordinate 
time framework, event horizon acts as an asymptotic line. 



\t 




FIG. 6: Radial motion of the test particle for the proper time, r (solid curve), and the coordinate time, t (dashed curve). These 
solutions are very similar to Schwarzschild case where in the proper time framework particles can cross the event horizon and 
in the coordinate time framework, event horizon acts as an asymptotic line. Both curves are plotted by using m = 0.05M, 
Q = 0.85M, q = O.OIM, £ = 23A/, E = 0.067M 



V. MOTION OF CHARGED PARTICLES ON THE REISSNER-NORDSTROM DE SITTER BLACK 

HOLE 

In a similar way than the previous case of the RNAdS black hole, we start with particles with angular momentum 
L 7^ 0, then radial quadrature is given by 

= [E~V^{r)][E-V{r)], (74) 
then we can write the basic equation of the trajectories as follows (A/3 = 1/A^ > 0) 



^2 



where the characteristic polynomial for the allowed trajectories is 



the coefficients are given by 



A2Q2i2 2A2ML2 ^2^i^ _ f^l 



^0 = — — , = — — — , CJ2 = -A l^ + Q - 



(75) 



6 

Pe{r)=J2uj,r^, (76) 
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Using the constant C 7^, we can write the angular quadrature as 

Now we will give a detailed description of trajectories as in the previous case. First we need to define the kind of 
motion in terms of his energy and angular momentum against the effective potential. 

• Orbit with angular momentum; here we have, for a specific combination of energy and angular momentum, the 
possibility to have planetary orbits, when the effective potential exhibits a minimum. Besides, we have three 
different kinds of orbits that depend on the energy of the particle, compared to the effective potential. The first 
kind orbits take place between the two maxima (that is planetary orbits). The second kind orbits have the same 
energy than those from the first kind, but they fall to the event horizon, r-|_. Third kind orbits are those with 
the same energy than the cases before, but they instead fall to the cosmologic horizon. r+-(_. 

• Radial Orbit (L = 0) this orbit corresponds to a free fall to cosmological horizon or event horizon. 




FIG. 7: Effective potential for charged particles on the geometry of RNdS for M = 1, Q = 0.85A/, £ = 400M, L = 0.33M^ 
q = 0.18M and m = 0.2M. 



A. Second Kind trajectory: Fall to event horizon 

In the de Sitter case there are trajectories that start from a finite distance, at rest, and fall to the event horizon. 
The motion is developed in the region tq > r > > po- From ()77|) with the negative sign and the characteristic 
polynomial Peir) = (i?o ^ r){re — r)^(ro — r)(r — po){r — do) we obtain the integral motion as follows, 

0(r) = C r =. (78) 

Jro {re - r)y/{Ro - r) (ro ~ r){r - po){r - ctq) 

Therefore the solution can be written as 

^0(r) = a, [PsF{^- ks) + 7sn((^; n,)] , (79) 
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where the constants are given by as = y Ps = -r^) ^"^^ Is = (ji'oi^°)(r°-ro) • ^l^o the parameters 

of the elhptic integrals are 



ks = \ 7^ w r, n, = — and Lp{r) = arcsm J- -— ^, (80) 

Fig. ([5]) shows us ahowed trajectories of the charged particles in this case. Here, we can observe that the trajectories 
starting from a distances bigger than event horizon are doomed to cross it and fall to singularity. 




FIG. 8: Second Kind motion: It corresponds to trajectory starting from a distances bigger than event horizon are doomed to 
cross it and fall to singularity. 



B. Third Kind Trajectory: Fall to cosmological horizon 



This case corresponds to the motion that starts at a finite distance from the event horizon and moves out from the 
field forces, crossing the cosmological horizon. In this case the allowed range is given by r > i?o- Integral Quadrature 
is 

r dr 
0W = C / =, (81) 

JRo {r - re)\/{r - i?o)(r - ro)(r - po){r - do) 

again the solution can be written in terms of elliptical Jacobi integrals 

^(b{r) = at[l3tF{ip;kt) + ^M^;h,nt)] , (82) 

where the constants are given by at = . ^ , Bt = ? — ^ — r and 7* = — . and the corresponding 

elliptic parameters are 



* \ [Ro- Po){ro - c^o) ' 



nt 



{Rq - <7o){re - rp) 
(''o - CTo)(re - Ro) ■ 



(84) 
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FIG. 9: Third Kind motion for a RNdS geometry: the particles moving out from the event horizon and then they are able to 
cross the cosmological horizon losing casual connection with physical observers. 



iPir) = arcsm J — ^, 85 

FIGlS] depicts these trajectories, that correspond to particles moving out from the event horizon and then they are 
able to cross the cosmological horizon losing casual connection with physical observers. 



C. Critical Trajectory 

As we have say, it is possible to obtain a class of orbits which posses two maxima at and Ri (r^ < Ri) with the 
same value of the constant E = En (see FIGlll). this case, the conditions 



V{r,) - V{R^) = Eij, (86) 

and 

are satisfied. In this scenario, we can identify three possible trajectories: particles falling from an initial distance, di 
(di < r_(-+), to the unstable circular orbit at Ri (first kind critical orbit); another one goes up to the unstable circular 
orbit at (second kind critical orbit) from an initial distance, d2 {d2 > r^); and finally, the particles going to an 
asymptotically circular orbits at Ri or (third kind critical orbit) depending on the sign of its initial velocity. Thus, 
the - (+) sign means that the particles falls (up) to the asymptotic orbit at {Ri). 



1. first kind critical orbit 



In this case, the particle starts its motion at di and falls to the asymptotic orbit at Ri. From ([TT)) with the negative 
sign and the characteristic polynomial given by Pg^r) = (r — i?o)^(r — ri)^(r — pi)(r — Ci), we obtain the solution 
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C 



^ In e(r, i?,)- In e(r,r,) 



an) 



where 



ri[u, V) 



(89) 



with 



r]{u,v) = u-v, and ^{u) = y/r]{u, Pi) r]{u,ai). 
In left panel of FIG. [10] we have shown a typical first kind trajectory. 



(90) 



2. second kind critical orbit 

This orbits have trajectories that go up asymptotically to the unstable circular orbit at from a distance d2 
(^2 > Thus, we must choose the plus sign in eq. (j77p with the characteristic polynomial given by Peif) = 

{Ri — r)'^(ri — r)'^(r — pi){r — ai), in which case we obtain the solution 



</>(2)(r) = 



C 



that we show in middle panel of FIG.fTUl 



In (-e(r, i?,)) - 77^ In (-^in n)) 



(91) 



3. third kind critical orbit 

Finally, when the particle is at an initial distance da, such that r.^ < da < it has the ability to go asymptotically 
to the unstable circular orbit at or Ri , depending on the sign of its initial velocity. Therefore, the solution can be 
write as 



0f(r)=± 



In (-e(r, i?,)) + 77^ In Q{r, n) 



^0, 



(92) 



where i^q is a constant of integration that fixes the initial conditions and the - (+) sign means that the particle 
approaches to the unstable circular orbit at rt (Ri). In right panel of FIG.[TU]we show this trajectories. 



VI. RADIAL MOTION 



The motion of radial particles {L = 0) is described by the equation 



drV 



qQ 



m^il + - 



(93) 



The effective potential in this case is showed in 

For different values of the constant of motion E, we have the following allowed orbits: The capture zone {E < Ec) 
and critical motion [E = E^). 
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FIG. 10: Critical trajectories for charged particles in RNdS spacetime. Left panel: First kind orbit. The particle starts its 
motion from di ((/)c^' = 0) and approaches asymptotically to unstable circular orbit at Ri; Middle panel: Second kind orbit. 
The particle starts its motion from d2 {(fi^^^ = 0) and approaches asymptotically to unstable circular orbit at r^. Right panel: 
Third kind orbit. The particle goes to the unstable circular orbit at n or Ri, depending on the sign of its initial velocity; All 
graphics has been plotted using M = 1, m = O.IM, Q = 0.85M, q = O.OIM, A = 400M 




FIG. 11: Effective potential for radial charged particles in the RNdS spacetime with m — 0.1, M = 1, Q = 0.8, q — 0.01, 
A = 400M 



A. Capture Zone 



In this case, we have that E < Ec, and from here, we study two kind of orbits: first kind orbit, which fall to the 
cosmological horizon, and second kind orbit, which fall to the event horizon. It is convenient to rewrite eq. (|93p in 
the form 



where the characteristic polynomial is 



dr\ _m-G{r)_^ ^^^^ 



dT A2 



G{r) =r'^ ^ ar^ +Pr-j, (95) 



and the constants are given by a = A^(l — E^/m'^), j3 — 2A^(Af — qQE/m?) and 7 = \^Q'^{1 — /ni^). 
Finally, we have the following relations in order to obtain the motion integrals 



± 



dr 



dT A 



(96) 
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(97) 



(98) 



First Kind: Fall to event horizon 



Radial motion shows more variety depending on the roots of the respective polynomial and the energy of the 
particles. This case corresponds to particles that start from a finite distance bigger than the last root of the effective 
potential and less than cosmological horizon. Here, allowed region corresponds to i?o > '''o > > c^i > > ^2, such 
that the function G{r) can be write as Gi{r) — [r — Ro){r — ro)(r — di){r — ^2), in which case, from eq. we 
obtain the following solution for the proper time 



where 



r(l)/'.^ - 



(99) 



ni = 



2Ro 



arcsm < 



{ro - r) 



\/{Ra - di){ro ~ ^2)' \'ni{Ro-r)' '^^ ^"^(^0-^2)' 

On the other hand, in the coordinate time framework the general solution corresponds to 

t^-^\r) = Aiti{r) + A2tii{r) + Astmir) + Aitiv{r), 
the divergent solution at the event horizon is tjj, which is given by 



{Ro 



and 



ni 



rp - di 
Ro — di 



(100) 



with 



tiiir) = {ro - r+)F{ipi;Ki) + {Rq - ro)n(?/'i; ki, ^2), 



(101) 



A2 = 



54 



1 



I , and 712 



' i?o - r+ ro - di 



(i?o - r+)(ro - r+)^(i?o - di){ro - (^2) ' ]j ro - r+ Ro - di 

FIG. [121 summarizes our result in this case and shows the fall to the event horizon. 



2. Second Kind: Fall to cosmological horizon 

Now, we are considering the radial motion of particles that fall to the cosmological horizon. In this case, the radial 
coordinate belong to the interval r++ > r > Rq. So that, the solution for the proper time is 



where 



T(2)(r) = r!2 



-^("02; Ki) + — — ^n(V'2; Ki, "-3) 



(102) 



no 



2rn 



^ \/{Ro- di){ro - d2)' 



. {r-Ro) , Ro-d2 

ip2 — arcsm < / — , and — 



n3{r - ro) ' 



''o - d2 
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In the coordinate time framework, after write the motion equation in terms of partial fractions, we obtain the 
general solution 



= A^tiir) + A2tiiir) + Astmir) + A^tivir), 
the divergent solution at the cosmological horizon is tj, which is given by 



with 



54 



ti{r) = (r++ - Ro)F{'>p2; Hi) + (Rq - ro)Il{^2; i^i^n^), 



and = 



(r++ - i?o)(r++ - rQ)^{Ro-di){ro-d2y 
FIG. [m shows the behavior of the radial fall to event horizon. 



r++ - Ro 



(103) 



(104) 



B. Radial Critical Trajectory 

One important case appears when a maximum in the effective potential exits. First at all, we consider that E < 1 
and 9 < 1, then we have four real roots, where the maximum occurs at r = Dc, 



where 



and the quadrature is 



Gii{r) = {r-Dcf9(r), 
T{r) = / '- — ^=dr. 



(105) 
(106) 

(107) 



Here, we can classify orbits in terms of the distance where the trajectory starts respect to Dc- We start with 
critical first kind, they correspond to trajectories that asymptotically move to Dc from greater distances i?o, and the 
region under consideration is given by i?o > ^ > Dc > Di > > D2. Then choosing negative sign for a fall from i?o, 
we obtain 



1 



(r - Dc) 



dr 



whose general solution is 



T{r) = ^[T/(r) +T//(r)], 



where the first integral is a regular function 

r/(r) = In 

second integral is divergent at Dc 



Dc + Rn + y/ffCRo) 
Dc + r+ yW) 



Ro-Dc 



r ~ Dc 



V2 



In 



Then, the final solution is given by 



A" Ac) 



Tf^(r) =ln 



Dc + Ro + VgiRo) 



Dc + r+ VsW 



1 , 

+ ^In 


'Ro-Dc 




V2 


r - Dc 


V2 



4Ro + 2V2v/.g(flo) 
4r + 2V2V5W 

4i?o + 2V2 V5(i?o) 



4r + 2\/2V5H 



(108) 



(109) 



(110) 



(111) 



(112) 
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The other kind corresponds to critical second kind trajectories where particles approximate to critical radios from 
distances less than Dq and take infinite time to reach that distance. 

If the particle starts at and we choose minus sign, the solution is given by 



A2 



.(<=) 



(r) = In 



Dc + ra + \Jg{rQ) 



1 



V2 



In 



Dc - r 



1 

71 



In 



4ro + 2V2^g{ra) 
Ar + 2V2\/5H 



(113) 




FIG. 12: Radial motion for charged particles in the coordinate time, t, and proper time, r with m = 0.1, M = 1, Q = 0.8, 
q = 0.01, A = 400M 



VII. FINAL REMARKS 



In this paper we studied in detail the motion of charged particles in the vicinity of the Reissner-Nordstrom (Anti)- 
dex Sitter black hole. We found, in the AdS case that exist different kinds of trajectories depending on the angular 
momentum of the particles. Specifically we found that exist 

• Planetary Orbit: In this case the orbit corresponds to a bounded trajectory that exhibit oscillation between 
two extremal distance (periastron and apastron). Besides, this trajectory have one particular solution called 
Reissner-Nordstrom limit. 

• Second Kind Trajectory: This trajectory is computed with the same parameter than planetary orbits. It 
corresponds to a trajectory that starts at rest from a finite distance less than periastron. This kind of motion 
represents the fall to the event horizon, and it is considered that have a turning point inside to cauchy horizon. 
Also it is shown an approximated solution called Reissner-Nordstrom limit 

• Critical Trajectories: There are trajectories of the first and second kind. The first one, starts at rest from a 
finite distance from the outside of the unstable circular orbit and then approximates asymptotically to it. The 
second kind approximates to the unstable orbit from a distance at the inner side of it. 

• Radial Trajectory: This trajectory has vanishing angular momentum and physically describes a radial fall to 
the event horizon starting at rest. 

• Circle Orbits: For some fixed values of the constants it is possible to find solutions of the equation of motion 
that represent stable and unstable circle orbits. 

In the dS case we mainly find two branches of classification that are summed as follow 

• Orbit with angular momentum: here we have for a specific combination of energy and angular momentum, the 
possibility to have planetary orbits, when the effective potential shows a minimum. Besides, we have three 
different kinds of orbits that depend on the energy of the particle compared with effective potential. First kind, 
corresponds to an orbit whose energy is less than the global maximum of the effective potential and starts at a 
large distance (bigger than the last root of the respective polynomial) and the particles are doomed to go to the 
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cosmological horizon. Second kind, are represented by orbits that start at a distance between the event horizon 
and the first root of the polynomial, they are doomed to fall to the event horizon. In this case the energy of 
the particle is also less than the global maximum of the effective potential. Third kind, corresponds to orbits 
whose energy is bigger than the effective potential and particles have possibilities to move to cosmological or 
event horizon depending on the initial conditions in the velocity. 

• Radial Orbit {I = 0) this orbit correspond to free fall to cosmological horizon or event horizon. Also, we can 
find the three kinds described for the orbit with angular momentum 

In sum we have characterized through a specialized study all the orbits for the black holes under consideration. We 
left the study of photons for a future work. 
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